ANALYSIS ARTICLE OPEN ACCESS 


DRUG DISCOVERY 


To Cite: 

Afolayan FID, Tarkaa CT. Network pharmacology-based 
assessment of anti-inflammatory action of phytocompounds 

erived from Nigella sativa and Moringa oleifera. Drug Discovery 2023; 
17: e13dd1016 

doi: https://doi.org/10.54905/disssi.v17i39.e13dd1016 


Author Affiliation: 
Cell Biology and Genetics, Department of Zoology, University of 


Ibadan, Ibadan, Nigeria 


“Corresponding Author 
Cell Biology and Genetics, Department of Zoology, University of 


Tbadan, Ibadan 
Nigeria 
Email: fidifede@gmail.com 


Peer-Review History 

Received: 03 January 2023 

Reviewed & Revised: 06/January/2023 to 03/March/2023 
Accepted: 07 March 2023 

Published: 12 March 2023 


Peer-Review Model 
External peer-review was done through double-blind method. 


Drug Discovery 
PISSN 2278-540X; elSSN 2278-5396 


© The Author(s) 2023. Open Access. This article is licensed under a 
Creative Commons Attribution License 4.0 (CC BY 4.0)., which permits use, 
sharing, adaptation, distribution and reproduction in any medium or 
format, as long as you give appropriate credit to the original author(s) and 
the source, provide a link to the Creative Commons license, and indicate if 
changes were made. To view a copy of this license, visit 
http://creativecommons.org/licenses/by/4.0/. 


« DISCOVERY 


\v .) SCIENTIFIC SOCIETY 


Drug Discovery 17, e13dd1016 (2023) 


Network pharmacology-based 
assessment of anti-inflammatory 
action of phytocompounds 
derived from Nigella sativa and 


Moringa oleifera 


Funmilayo ID Afolayan’, Christopher T Tarkaa 


ABSTRACT 


Studies have shown that Nigella sativa and Moringa oleifera possess anti- 
inflammatory activity, however, the molecular targets involved in their 
mechanisms of action are not known. A network-based pharmacology analysis 
was done to predict molecular targets and phytocompounds involved in the anti- 
inflammatory activity of N. sativa and M. oleifera. Phytocompounds of N. sativa 
and M. oleifera were retrieved from Dr Duke's Phytochemical and Ethnobotanical 
databases and Indian Medicinal Plants, Phytochemistry and Therapeutics 
database. Target proteins were obtained from Binding DB. A compound-target- 
pathway network was constructed with the Cytoscape tool and a network of 
protein-protein interactions was established with the STRING database. Lead 
proteins identified from the compound-target-pathway network were further 
studied for their interactions with N. sativa and M. oleifera by molecular docking. 
Similarly, biological pathways involved in the anti-inflammatory activity of the 
phytocompounds were identified with the Kyoto Encyclopedia of Genes and 
Genomes (KEGG) analysis. Forty-two and twenty-seven bioactive compounds 
from N. sativa and M. oleifera respectively were successfully identified 
corresponding to 98 targets that were screened out for anti-inflammation. Based 
on network pharmacology data and molecular docking, beta-sitosterone, 
astragalin, beta-amyrin and quercetin were predicted to inhibit key inflammatory 
proteins. Target proteins, including NR1H3, F2, AKT1, HSP90AA1, IL2, NFKB1, 
PTGS2, ALOX5, TNF-a, IL6 and IFN-y, as well as signaling pathways such as 
TNF, MAPK, IL-17 and HIF-1 were linked to the anti-inflammatory activity of N. 
sativa and Moringa oleifera. Functional enrichment analysis predicted that 33 
inflammatory pathways were modulated by N. sativa and M. oleifera which need 


further analyses for confirmation. 
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1. INTRODUCTION 


Inflammation is the body’s defense mechanism that eliminates injurious stimuli such as irritants, damaged cells, infection and other 
stimulants and initiates the process of healing (Jain et al., 2015). Acute inflammation, with exudation of fluid and plasma proteins at 
the affected site, intravascular activation of platelets and the emigration leukocytes as its main characteristics occurs very fast and 
the process can last for a few minutes to numerous days. Chronic inflammation on the other hand occurs when acute inflammation 
occurs repeatedly or continuously, with the process lasting for several weeks to months and even years (Kosala et al., 2018; Kumar 
et al., 2018). Inflammation is like a two-edged sword because although it removes injurious stimuli and initiates the healing process, 
dysregulated inflammation can lead to tissue and cell injury, chronic infection, chronic diseases and cancers (Roth, 2012). Some 
drugs are presently used to manage inflammatory conditions. These drugs include steroidal anti-inflammatory drugs (SAIDs) and 
non-steroidal anti-inflammatory drugs (NSAIDs). Non-steroidal anti-inflammatory drugs remain a mainstay of therapy for the 
treatment of a large number of inflammatory diseases. Non-steroidal anti-inflammatory drugs do not cure the disease but instead 
modify the inflammatory response to the disease. The use of NSAIDs is frequently limited by gastrointestinal side effects, that range 
from dyspepsia to bleeding from ulceration which can be life-threatening. Non-steroidal anti-inflammatory drugs have been 
reported to inhibit the cyclooxygenase (COX) pathway (Borquaye et al., 2020). 

Cyclooxygenase is required for the conversion of arachidonic acid to prostaglandins, prostacyclin and thromboxanes. 
Specifically, prostaglandin causes vasodilation, increases the temperature set-point in the hypothalamus and play a role in anti- 
nociception, prostacyclins cause vasodilation, inhibition of platelet aggregation, relaxation of gastrointestinal and uterine muscle 
and release of pituitary hormones and thromboxane plays a role in platelet adhesion (Jain et al., 2015; Ghlichloo and Gerriets, 2021). 
Two cyclooxygenase isoenzymes, cyclooxygenase-1 (COX-1) and cyclooxygenase-2 (COX-2) are known. Cyclooxygenase-1 is 
constitutively expressed in the body, and it plays an important role in the maintenance of gastrointestinal mucosa lining, kidney 
function and platelet aggregation. Cyclooxygenase-2 is not constitutively expressed in the body, rather, it is expressed during an 
inflammatory response. The Majority of the NSAIDs are nonselective and are known to inhibit both COX-1 and COX-2. A few 
NSAIDs such as celecoxib and refecoxib are selective and can target only COX-2, thus having a different side effect profile 
(Ghlichloo and Gerriets, 2021). Other side effects of anti-inflammatory drugs include gastric ulcers, gastric irritation, alterations in 
renal function, hepatic injury, effects on blood pressure and platelet inhibition which may result in increased bleeding (Geremew et 
al., 2015). It is therefore imperative to develop newer anti-inflammatory drugs that are safer for use by humans. 

Traditional medicines have been used for ages to maintain health and improve various diseases including chronic inflammatory 
diseases. These agents comprise various plant sources and are associated with minimal side effects (compared to modern allopathic 
medicines), show synergistic effects by interacting with multiple proteins and modulating multiple pathways (Yuan et al., 2017). 
About 25% of medicines used today are developed from traditionally used medicines; however, there is a paucity of scientific 
evidence for the mechanism of multicomponent treatments (Duyu et al., 2020). Nigella sativa is a plant that belongs to the family 
Ranunculaceae and it is an annual herbaceous plant that is native to Southern Europe, Southwest Asia and North Africa. The plant 
is commonly known as black cumin and it is cultivated in many countries including South Europe, India, Syria, Turkey, Pakistan, 
the Middle Eastern Mediterranean region and Saudi Arabia (Hannan et al., 2021). N. sativa has been studied extensively for its 
biological activities and therapeutic potential and shown to possess a wide spectrum of activities such as diuretic, antihypertensive, 
anticancer, anti-inflammatory, antidiabetic, immunomodulatory, anthelmintic, antimicrobial, hepatoprotective, renal protective, 
gastroprotective, spasmolytic, bronchodilator and antioxidant activities (Kooti et al., 2016). 

Moringa oleifera is a plant that belongs to the family Moringaceae and it is a small tree of the sub-Himalayan region of North 
West India, that is now indigenous to many regions in Africa, South East Asia, Arabia, the Caribbean and Pacific Islands and South 
America. The plant is commonly known as moringa and it is widely used for its health benefits. Among commoners, it has earned 
its name as the miracle tree because of its amazing healing abilities for various ailments (Razis et al., 2014). Each part of the plant 
(fruits, seeds, leaves, flowers, bark and roots) has one or more benefits and are being used traditionally for different purposes, 
making it one of the most used plants in the world today (Mallenakuppe et al., 2019). M. oleifera has also been studied extensively 
for its biological activities and therapeutic potential and has also shown to possess a wide spectrum of activities such as 
antidiabetic, anticonvulsant, antidepressant, anti-inflammatory, anticancer, immunomodulatory, antihyperlipidemic, antimicrobial, 
antiulcer, antiarthritic, antioxidant, SCAR marker, antimutagenic, antifungal, antiurolithiatic activities (Srivastava et al., 2020). 

Network pharmacology is a research field that uses several components and several targets as the starting point, to study the 
molecular mechanism of action of various drugs and plants in the treatment of diseases (Zhang et al., 2020). This method integrates 
system biology, multiple pharmacology and computer analysis technology. Unlike the previous concept of single-component/single 
target, network pharmacology uses a large number of databases and statistical algorithms to identify the synergistic effects of 
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several components, several targets and several pathways of disease (Yu et al., 2020). It also explores the interaction between the 
drugs and potential targets and establishes a “compounds-targets-pathways” network to associate drugs systematically and 
comprehensively (Xie et al., 2015). Network pharmacology is thus a more holistic approach in carrying out research than the single- 
component/single-target approach. 

Molecular docking is usually used to validate the interactions between the compounds and target (Que et al., 2021). Molecular 
docking is a high-throughput technology that is widely applied to study the active sites of drugs by simulating interactions between 
receptors and drug molecules (Yu et al., 2020). It has therefore played an immense role in the research of natural products (Gogoi et 
al., 2017). In this study, network pharmacology was used to construct the active components-targets-pathways network diagram of 
phytocompounds present in Nigella sativa and Moringa oleifera. Additionally, molecular docking technology was used to predict the 
potential active components and mechanism of action of the plants. These findings set the basis for further research in this direction. 


2. MATERIALS AND METHODS 


Collection of Active Components 

The reported phytoconstituents from Nigella sativa and Moringa oleifera were retrieved from Dr Duke’s Phytochemical and 
Ethnobotanical Database (DPED) (https://phytochem.nal.usda.gov/phytochem/search) and Indian Medical Plants, Phytochemistry 
and Therapeutics (IMPPAT) (https://cb.imsc.res.in/imppat/home) using the keywords “Nigella sativa” and “Moringa oleifera” 
respectively (Gu et al., 2013). Canonical SMILES of each compound along with their molecular weight and the molecular formula 
were recorded from PubChem (https://pubchem.ncbi.nlm.nih.gov/) and Chemspider (https://www.chemspider.com/) databases 
(Kim et al., 2015). 


Prediction and Screening of the Inflammatory Targets of the Active Components 

The targets of each compound were queried using Binding DB (https://www.bindingdb.org/) at a percentage similarity of 70%. The 
gene code of each target was identified using UniProt (https://www.uniprot.org/). Similarly, targets for inflammation were 
identified with reference to the Gene Cards database (https://www.genecards.org/) (Zhang et al., 2020). 


Drug likeness and ADMET Profile 

The drug-likeness score for each compound targeting proteins related to inflammation was predicted using the Mol Soft database 
(https://molsoft.com/) based on “Lipinski’s Rule of Five” model (Duyu et al., 2020). Likewise, the probability for human intestinal 
absorption, caco-2 and blood-brain barrier permeability, human oral bioavailability, P-glycogen substrate, P-glycoprotein, CYP3A4, 
CYP2C9, CYP2D6 and CYP1A2 inhibition, eye irritation, Ames mutagenesis, human ethera-go-go inhibition, carcinogenicity, 
hepatotoxicity, etc. were predicted using admetSAR2.0 (http://Immd.ecust.edu.cn/admetsar2/) (Yang et al., 2019). 


Pathway Enrichment Analysis 

The potential gene targets were uploaded to the KOBAS 3.0 database (http://kobas.cbi.pku.edu.cn/kobas3) for KEGG pathway 
annotation. These enriched pathways were considered significant ones that will be affected by Nigella sativa and Moringa oleifera in 
the treatment of inflammation (Yu et al., 2020). 


Construction of Protein-Protein Interaction Network 

The gene targets were imported into the STRING database (https://string-db.org/) to generate the target protein interaction 
network. The TSV file was downloaded and imported into Cytoscape 3.8.2 software and the network was analyzed using the 
Network Analysis function (Ye et al., 2020). 


Construction of components-targets—pathways network 

The active components, gene targets and enrichment pathways were imported into the Cytoscape 3.8.2 software to construct the 
components-targets-pathways network. Active components, gene targets and enrichment pathways represent the input nodes. The 
association between the 2 nodes is represented by an edge (Yu et al., 2020). 


Molecular Docking 
Based on the network pharmacological outcome, molecular docking analyses were carried out between the phytoconstituents and 


AKT1 and HSP90AA1 (the proteins that interact with most targets) and NR1H3 and F2 (the proteins that are targeted by most 
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phytoconstituents). Because IL-2, NFKB1, PTGS2 (COX-2), ALOX5, TNF, IL-6 and IFN-y are closely linked to inflammation, they 
were also subjected to docking studies (Yu et al., 2020). The 3D structures of the proteins were obtained from the protein data bank 
(https://www.rcsb.org/). The PDB ID of AKT1, HSP90AA1, NR1H3, F2, IL-2, NFKB1, PTGS2, ALOX5, TNF-a and IL-6 are 1UNR, 
3OWB, 2ACL, 1PPB, 1M47, 1A3Q, 5IKT, 308Y, 1TNF, 1ALU and 1FG9 respectively. To avoid docking interference, heteroatoms 
and water molecules were removed with discovery studio 2021 (Duyu et al., 2020) and hydrogen bonds were added using Chimera 
1.14 Johnson et al., 2020). The active sites of the proteins were obtained from the CASTp 3.0 database (http://sts.bioe.uic.edu/castp/) 
(Tian et al., 2018). The SDF (structure-data file) structures of the phytoconstituents and standard ligand i.e., diclofenac was obtained 
from the PubChem database. The ligands were docked to the protein targets at their active sites and the binding affinities were 
determined using Auto dock Vina from PyRx. Molecular interactions between proteins and ligands were viewed with Chimera 1.14 
and discovery studio 2021 (Johnson et al., 2020). 


3; RESULIS 


Active components and their targets 

Eighty-one (81) phytoconstituents were identified from N. sativa, out of which forty-two (42) were associated with inflammation 
(Table 1). Also, fifty-seven (57) were identified from M. oleifera in which twenty-seven (27) modulated proteins involved in 
inflammation (Table 2). Diclofenac was chosen as the standard drug because it is one of the common non-steroidal anti- 


inflammatory drugs used in the treatment of inflammation and pain (Kolodziejska and Kolodziejczyk, 2018). 


Drug-likeness and ADMET profile 

Out of Forty-two phytoconstituents from N. sativa, twenty-one phytoconstituents were predicted to score positive and negative 
drug-likeness hit each. Similarly, out of twenty-seven compounds from M. oleifera, seventeen and ten were predicted to score 
positive and negative drug likeness hit respectively. Rutin was predicted to score the highest drug-likeness score of 0.91 for N. sativa 
while delta-5-avenasterol was predicted to score the highest drug likeness score of 0.85 for M. oleifera. 

Twenty-four phytoconstituents from N. sativa were predicted to violate 1 Lipinski rule, two phytoconstituents (Caspase 8 
inhibitor and Quercetin-3’-glucoside) were predicted to violate 2 Lipinski rules and one phytoconstituent (Rutin) was predicted to 
violate three Lipinski rules. Fourteen phytoconstituents from M. oleifera were predicted to violate one Lipinski rule and one 
compound (Maltotriose) was predicted violate three Lipinski rules. 

The list of phytoconstituents and drug likeness scores are summarized in tables 3 and 4. The ADMET profile of each 
phytoconstituent is represented in the heat maps (Figures 1 and 2). All the phytoconstituents from N. sativa are predicted to have 
positive human intestinal absorption scores. All the phytoconstituents from N. sativa except Astragalin, Betulic acid, Hederagenin, 
Quercetin-3’-glucoside and Rutin are blood-brain barrier permeants, twenty-two phytoconstituents have positive bioavailability 
scores and all phytoconstituents apart from Caspase 8 inhibitor and Carvacrol are predicted to be non-inhibitors of P-glycogen. 
Caspase 8 inhibitor, Dithymoquinone and Nigellone are predicted to be inhibitors of CYP3A4 and nine phytoconstituents are 
predicted to be inhibitors of CYP1A2. Twelve phytoconstituents are predicted to cause eye irritation. Astragalin, Butyrospermol, 
Quercetin-3’-glucoside and Rutin are predicted to be mutagenic, eight phytoconstituents are predicted to be inhibitors of human 
ether a-go-go, eight phytoconstituents showed a tendency towards hepatotoxicity and none of the phytoconstituents showed a 
tendency towards carcinogenicity. 


Table 1 Phytoconstituents of Nigella sativa and the inflammatory proteins they target 


S/N_ | Phytoconstituents Targets 

PRKAB1, CACNAIC, F7, F10, GRIN2B, GRIN1, ITGB3, PTPRC, 
ELANE, NFKB1, RORC, NR1H3, NR1H2, PTPN1, F2 

2 24-methylene-cycloartanol | NR1H3 

PRKAB1, CACNAIC, F7, F10, GRIN2B, GRIN1, ITGB3, PTPRC, 
ELANE, NFKB1, RORC, NR1H3, NR1H2, PTPN1, F2 

CYP24A1, GRIN2B, GRIN1, ITGB3, RORC, NR1H3, NR1H2, F2, 
RXRA 
PRKAB1, CACNAIC, F10, GRIN2B, GRIN1, ITGB3 
RORC, NR1H3, NR1H2, PTPN1, F2 

6 Astragalin ALOXS5, IL2, MMP9, NOX4, PTPN1, F2, RPS6KA3, TNF 


1 24-ethyl-lophenol 


3 24-methyl-lophenol 


4 5-dehydro-avenasterol 


5 Alpha-spinasterol 
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. F7, F10, GRIN2B, GRIN1, ITGB3, PTPRC, ELANE, NFKB1, 
7 Beta-amyrin 
RORC, NR1H3, NR1H2, PTPN1, F2 
8 Betulic acid F7 
PRKAB1, CACNAIC, F7, F10, GRIN2B, GRIN1, ITGB3, PTPRC, 
9 Butyrospermol 
ELANE, NFKB1, RORC, NR1H3, NR1H2, PTPN1, F2 
PRKABI1, F10, GRIN2B, GRIN1, ITGB3, RORC, NR1H3, NR1H2, 
10 Campesterol 
PTPNI, F2 
11 Carvacrol ATP2A1 
12 Carvone BIRC2 
ACE, CAPN2, CAPN1, CASP3, CASP7, CASP8, CASP9, CTSK, 
13 Caspase 8 Inhibitor CTSS, ITGA2B, ITGB3, ELANE, MALT1, PSEN1, CTSL, REN, 
CMA1 
CAPN2, CAPN1, CASP3, CASP7, CASP8, CASP9, 
14 Caspase inhibitor VI CTSK, CTSS, NFKBIA, ITGA2B, ITGB3, ELANE, MALT1, 
PSEN1, CTSL 
PRKABI1, F10, GRIN2B, GRIN1, ITGB3, RORC, NR1H3, NR1H2, 
15 Cholesterol 
PTPNI, F2 
PRKAB1, CACNAIC, F7, F10, GRIN2B, GRIN1, ITGB3, PTPRC, 
16 Citrostadienol 
ELANE, NFKB1, RORC, NR1H3, NR1H2, PTPN1, F2 
17 Cycloeucalenol NR1H3 
18 Damascenine ALOX5 
19 Dehydroascorbic-acid ALOX5 
PRKAB1, CACNAIC, F10, GRIN2B, GRIN1, ITGB3 
20 Delta 7-avenasterol 
RORC, NR1H3, NR1H2, PTPN1, F2 
21 Dithymoquinone BIRC2 
22 Dl-alanine-15n GRIN1 
23 Eicosadienoic-acid PTPN7, PPARA 
PRKAB1, CACNAIC, F7, F10, GRIN2B, GRIN1, ITGB3, PTPRC, 
24 Gramisterol 
ELANE, NFKB1, RORC, NR1H3, NR1H2, PTPN1, F2 
. F7, F10, PTPRC, ELANE, NOS2, NFKB1, RORC, NR1H3, 
25 Hederagenin 
PPARA, PTGS2, PTPN1 
26 Indole-3-acetic-acid C5, PTGS2 
PRKAB1, CACNAIC, F7, F10, GRIN2B, GRIN1, ITGB3, PTPRC, 
27 Lophenol 
ELANE, NFKB1, RORC, NR1H3, NR1H2, PTPN1, F2 
28 N-acetyl-L-cysteine ACE 
29 Nigellone BIRC2 
oe CACNAIC, F10, GRIN2B, GRIN1, ITGB3, RORC, NR1H3, 
30 Obtusifoliol 
NR1H2, PTPN1, F2 
31 Octadeca-9,12-dienoic acid | PTPN7, PPARA 
32 Oleic acid PTPN7, PPARA 
PRKAB1, F10, GRIN2B, GRIN1, ITGB3, RORC, NR1H3, NR1H2, 
33 Phytosterols 
PTPNI, F2 
. : ALOX5, IL2, MMP9, NOX4, PTPN1, F2, RPS6KA3, TXNRD1, 
34 Quercetin-3'-glucoside 
TNF 
ALOX5, IL2, MMP9, NOX4, PTPN1, F2, RPS6KA3, TXNRD1, 
35 Rutin 
TNF 
PRKAB1, CACNAIC, F10, GRIN2B, GRIN1, ITGB3, RORC, 
36 Stigmast-7-en-3-beta-ol 
NR1H3, NR1H2, PTPN1, F2 
37 Stigmastanol RORA, NR1H3 
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F7, F10, GRIN2B, GRIN1, ITGB3, PTPRC, ELANE, NFKB1, 
38 Taraxerol 
RORC, NR1H3, NR1H2, PTPNI, F2 
39 Thymohydroquinone ATP2A1 
40 Thymol ALOX5, ATP2A1 
41 Thymoquinone ALOX5 
CACNAIC, F10, GRIN2B, GRIN1, ITGB3, RORC, NR1H3, 
42 Tirucallol 
NR1H2, PTPN1, F2 
43 *Diclofenac CXCRI1, PTGS2, CXCL8 
*Standard drug 


Table 2 Phytoconstituents of Moringa oleifera and the Inflammatory proteins they target 


S/N | Phytoconstituents Targets 
1 24-Methylene-cholesterol | CACNAIC, HSP90AA1, BIRC2, NR1H3 
2 28-Isoavenasterol acetate RORC, F10, PTPN1, NR1H2, NR1H3 
3 4-hydroxymellein IL6, PRKCA 
4 Benzyl-amine F2, F10 
5 Beta-sitosterone CACNAIC, HSP9IDAAI, BIRC2, NR1H3 
; ITGB3, F2, GRIN2B, PRKAB1, GRIN1, RORC, F10, PTPN1, 
6 Brassicasterol 
NR1H2, NR1H3 
7 Caffeic-acid LCK, MMP2, MMP1, MMP9, EGFR 
Campestanol RORA, NR1H3 
PRKABI1, F10, GRIN2B, GRIN1, ITGB3, RORC, NR1H3, NR1H2, 
9 Campesterol 
PTPNI1, F2 
PRKAB1, F10, GRIN2B, GRIN1, ITGB3, RORC, NR1H3, NR1H2, 
10 Clerosterol 
PTPNI, F2 
PRKAB1, F10, GRIN2B, GRIN1, ITGB3, RORC, NR1H3, NR1H2, 
11 Delta-5-avenasterol 
PTPNI1, F2 
PRKAB1, F10, GRIN2B, GRIN1, ITGB3, RORC, NR1H3, NR1H2, 
12 Delta-7-avenasterol 
PTPNI1, F2, CACNAIC 
PTGS2, NDUFA4, AKT1, ALOX5, PPARA, GSTP1, 
13 | Delta-tocopherol 
PTGER3 
HDAC1, ABL1, HDAC10, VEGFA, FGF2, ABL2, HDAC3, HDAC2, 
14 | D-galactose 
CDK2, HDACI, HDAC6, HDAC4, FGF1, CDK1, HDAC8 
15 Dl-alanine-15n GRIN1 
16 | Dl-arginine NOS3, NOS2 
‘a ae HDAC1, ABL1, HDAC10, VEGFA, FGF2, ABL2, HDAC3, HDAC2, 
ea HDAC9, HDAC6, HDAC4, HDAC5, FGF1, HDAC8 
18 Gamma-tocopherol PTGS2, AKT1, PTGER4, ALOX5, PPARA, GSTP1, PTGER3 
19 Indole-3-acetic-acid C5, PTGS2 
20 Indoleacetonitrile MPO, AHR 
PGF, PTGS2, IGFIR, BACE1, F2, LCK, MET, FLT3, VEGFA, TERT, 
GSK3B, ALOX5, EGFR, PIK3CG, CDK6, PPARA, CSNK2A1, 
21 Kaempferol 
YP1B1, PTPN1, NOX4, AHR, CDK2, NTRK2, CDK5, PIM1, CDK1, 
KIT, ABCC1 
HDAC1, ABL1, HDAC10, VEGFA, FGF2, ABL2, HDAC3, HDAC2, 
22 L-rhamnose 
CDK2, HDACI, HDAC6, HDAC4, FGF1, CDK1, HDAC8 
HDAC1, PRKCA, ABL1, HDAC10, VEGFA, FGF2, ABL2, HDAC3, 
23 Maltotriose HDAC2, CDK2, HDACI, HDAC6, HDAC4, HDAC5, FGF1, CDK1, 
HDAC8 
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24 | Phytosterols 


PRKABI1, F10, GRIN2B, GRIN1, ITGB3, RORC, NR1H3, NR1H2, 
PTPNI, F2 


25 Quercetin 


PGF, PTGS2, IGFIR, BACE1, F2, LCK, MET, FLT3, VEGFA, TERT, 
GSK3B, ALOX5, EGFR, PIK3CG, CDK6, PPARA, CSNK2A1, 
YP1B1, PTPN1, NOX4, AHR, CDK2, NTRK2, CDK5, PIM1, CDK1, 


KIT, ABCC1 
26 Stigmastanol RORA, NR1H3 
27 Vanillin ERN1 
28 *Diclofenac CXCR1, PTGS2, CXCL8 


*Standard drug 


Table 3 The molecular formula, physicochemical properties and drug-likeness scores of phytoconstituents from N. sativa 


Molecular 
Compounds Molecular formula NHBA | NHBD | MolLogP | DLS 
mass (g/mol) 
24-ethyl-lophenol C30Hs20 428.40 1 1 8.24 0.40 
24-methylene-cycloartanol | C3i1Hs20 440.40 1 1 8.49 -0.48 
24-methyl-lophenol C20H500 414.39 1 1 7.79 0.17 
5-dehydro-avenasterol CooHasO 410.35 1 1 7.34 0.57 
Alpha-spinasterol C20H4sO 412.37 1 1 7.87 -0.05 
Astragalin CaH2001 448.10 11 7 -0.12 0.67 
Beta-amyrin C30H500 426.39 1 1 7.95 -0.22 
Betulic acid C30H4sOs 456.36 3 2 7.05 0.25 
Butyrospermol C30H500 426.39 1 1 8.81 0.46 
Campesterol C2sH4sO 400.37 1 1 7.87 0.59 
Carvacrol CioHuO 150.10 a 1 3.44 -0.35 
Carvone CioHuO 150.1 1 0 2.99 -1.40 
Caspase 8 inhibitor C30Ha3FN40O11 654.29 11 5 1.72 -1.33 
Caspase inhibitor VI C21H2sFN3O7 453.19 7 4 1.02 -1.58 
Cholesterol C27H460 386.35 1 1 7.44 0.49 
Citrostadienol C30Hs00 426.39 1 1 7.74 0.47 
Cycloeucalenol C30H500 426.39 1 1 7.22 0.25 
Diamascenine CioHi3NOs 195.09 3 1 1.36 -0.86 
Dehydroascorbic-acid CoH6Oc 174.02 6 2 -1.95 -0.60 
Delta 7-avenasterol C20oHasO 412.37 1 1 7.95 0.25 
Dithymoquinone CooH24O4 328.17 4 0 2.76 -0.89 
Dl-alanine-15N C3H7NO2 89.05 3 3 -2.64 -0.96 
Eicosadienoic-acid C20H36O2 308.27 2 1 7.61 -0.30 
Gramisterol C20H4sO 412.37 a 1 7.39 0.13 
Hederagenin C30HasO4 472.36 4 3 4.99 0.47 
Indole-3-acetic-acid CioH9NO2 175.06 2 2 1.62 -1.43 
Lophenol CosH4sO 400.37 1 1 7.36 0.11 
N-acetyl-L-cysteine CsH»9NOsS 163.03 + 3 -0.92 0.30 
Nigellone C20H2404 328.17 4 0 2.76 -0.89 
Obtusifoliol C30H500 426.39 1 1 7.38 0.85 
Octadeca-9,12-dienoic acid | CisH32O2 280.24 2 1 6.60 -0.30 
Oleic acid CisH3402 282.26 2 1 7.11 -0.30 
Phytosterols C2H500 414.39 1 1 8.45 0.78 
Quercetin-3'-glucoside C21 H20012 464.10 12 8 -0.54 0.68 
Rutin Co7H300 16 610.15 16 10 -1.55 0.91 
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Stigmast-7-en-3-beta-ol C20H500 414.39 1 1 8.58 0.18 
Stigmastanol C2H320 416.40 ai 1 8.83 0.22 
Taraxerol C30Hs00 426.39 1 1 8.11 -0.90 
Thymohydroquinone CioH1sO2 166.10 2 2 2.34 -1.28 
Thymol CioHuO 150.10 1 1 3.43 -0.54 
Thymoquinone Ci0H1202 164.08 2 0 2.07 -1.03 
Tirucallol C30Hs00 426.39 1 1 8.74 0.55 
*Diclofenac CiusHuCl2aNO2 295.02 2 2 4.52 0.59 
NHBA- Number of hydrogen bond acceptors, NHBD- Number of hydrogen bond donors, DLS- Drug-likeness score 
Table 4 The molecular formula, physicochemical properties and drug-likeness scores of phytoconstituents from M. oleifera 
Compounds Molecdériomats | | eeea. | NED: | meitoge [OES 
mass (g/mol) 
24-methylene-cholesterol | C2sH1cO 398.35 1 1 7.47 0.51 
28-isoavenasterol acetate | Cs2Hs202 468.40 2. 0 8.25 1.20 
4-hydroxymellein Ci0H1004 194.06 4 2 1.31 -0.01 
Benzyl-amine C7HoN 107.07 1 2 1.04 -1.33 
Beta-sitosterone C2oHasO 412.37 1 0 8.40 0.68 
Brassicasterol CosHasO 398.35 1 1 7.08 0.34 
Caffeic-acid CoHsO« 180.08 4 3 1.27 -0.35 
Campestanol C28sH500 402.39 1 1 8.25 0.02 
Campesterol C2sHasO 400.37 1 1 7.87 0.59 
Clerosterol C29H4s0 412.37 1 1 8.58 0.68 
Delta-5-avenasterol CooH4sO 412.37 1 1 7.82 0.85 
Delta-7-avenasterol C29H4s0 412.37 1 1 7.95 0.25 
Delta-tocopherol Co7HasO2 402.35 2 1 9.20 0.33 
D-galactose CoH1206 180.06 6 5 -3.02 -0.12 
Dl-alanine-15n CsH7NO2 89.05 3 3 -2.64 -0.96 
Dl-arginine CoH1aNaO2 174.11 4 7 -3.93 -0.16 
D-xylose CsH100s 150.05 5 4 -2.93 -1.37 
Gamma-tocopherol CosHasO2 416.37 2 1 9.48 0.48 
Indole-3-acetic-acid CioHoN O2 175.06 2 2 1.62 0.50 
Indoleacetonitrile CioHsN2 156.07 1 1 1.86 -2.06 
Kaempferol C15H1006 286.05 6 4 1.61 0.11 
L-rhamnose CoH120s 164.07 5 4 -2.30 -1.05 
Maltotriose CisH32016 504.17 16 11 -5.06 0.11 
Phytosterols C2H500 414.39 1 1 8.45 0.78 
Quercetin CisH1007 302.04 7 5 1.19 0.52 
Stigmastanol C2Hs20 416.40 1 1 8.83 0.22 
Vanillin CsHsO3 152.05 3 1 0.96 -1.24 
*Diclofenac CuHuCl2NO2 295.02 2 2 4.52 0.59 
NHBA- Number of hydrogen bond acceptors, NHBD- Number of hydrogen bond donors, DLS- Drug-likeness score 


Drug Discovery 17, e13dd1016 (2023) 8 of 31 


ANALYSIS ARTICLE OPEN ACCESS 


Human intestinal absorotion 
P-glycogen inhibitor 

CYP2C9 inhibition 

CYP2D6 inhibition 

CYP1A2 inhibition 

Eye irritation 
Carcinogenicity 


ir mutagenesis 


24-ETHYL-LOPHENOL 
24-METHYLENE-CYCLOARTANOL 
24-METHYL-LOPHENOL 
5-DEHYDRO-AVENASTEROL 
ALPHA-SPINASTEROL 
ASTRAGALIN 

BETA-AMYRIN 

BETULIC ACID 

BUTYROSPERMOL 
CAMPESTEROL 

CARVACROL 

CARVONE 

CASPASE 8 INHIBITOR 

CASPASE INHIBITOR VI 
CHOLESTEROL 
CITROSTADIENOL 
CYCLOEUCALENOL 
DAMASCENINE 
DEHYDROASCORBIC-ACID 


a 
DELTA 7-AVENASTEROL as | li 
DITHYMOQUINONE EE me 
DL-ALANINE-15N Ei 
EICOSADIENOIC-ACID 
GRAMISTEROL 
HEDERAGENIN 
INDOLE-3-ACETIC-ACID 
LOPHENOL 
N-ACETYL-L-CYSTEINE 
NIGELLONE 
OBTUSIFOLIOL 
OCTADECA-9,12-DIENOIC ACID 
OLEIC ACID 
PHYTOSTEROLS 
QUERCETIN-3'-GLUCOSIDE 
RUTIN 
STIGMAST-7-EN-3-BETA-OL 
STIGMASTANOL 
TARAXEROL 
THYMOHYDROQUINONE 
THYMOL 
THYMOQUINONE 
TIRUCALLOL 
DICLOFENAC 


50 


Percentile 


Figure 1 Heat map of ADMET profile of phytoconstituents of Nigella sativa 
Green, yellow and red boxes indicate good, average and bad phytoconstituents for a particular parameter. 
24-ethyl-lophenol, 24-methyl-lophenol, cholesterol, lophenol, N-acetyl-4-cysteine, phytosterols, stigmat-7-en-beta-ol and stigmastanol had perfect positive outcomes. 


All the phytoconstituents from M. oleifera apart from D-galactose, DL-arginine, D-xylose, L-rhamnose and Maltotriose are 
predicted to have positive human intestinal absorption scores. All the phytoconstituents apart from 4-hydroxymellein, Caffeic acid, 
Kaempferol and Quercetin are blood-brain barrier permeants, eighteen phytoconstituents have positive bioavailability scores and 
all the compounds apart from 28-isoavenasterol acetate, Beta-sitosterone, Caffeic acid and Delta-5-avenasterol are predicted to be 
non-inhibitors of P-glycogen. Kaempferol and Quercetin are predicted to be inhibitors of CYP3A4, Benzyl-amine, Indoleacetonitrile 
and Kaempferol are predicted to be inhibitors of CYP2C9, Indoleacetonitrile is predicted to be an inhibitor of CYP2D6, 4- 
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hydroxymellein, Benzyl-amine, Indoleacetonitrile, Kaempferol and Quercetin are predicted to be inhibitors of CYP1A2. Nine 


phytoconstituents are predicted to cause eye irritation. 4-hydroxymellein, DL-arginine, D-xylose, Kaemferol and Quercetin are 
predicted to be mutagenic, eight phytoconstituents are predicted to be inhibitors of human ether a-go-go, Indole-3-acetic acid, 
Indoleacetonitrile, Kaempferol and Quercetin showed a tendency towards hepatotoxicity and Benzyl-amine showed tendency 


towards carcinogenicity. 
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Figure 2 Heat map of ADMET profile of phytoconstituents of Moringa oleifera 


Green, yellow and red boxes indicate good, average and bad phytoconstituents for a particular parameter. 


Campestanol and stigmastanol had perfect positive outcomes. 


Enrichment/network analysis 
Gene set enrichment analysis identified thirty-three (33) modulated inflammatory pathways. The predicted pathways in 
inflammation having the highest count of the gene set were MAPK signaling pathway and cytomegalovirus infection with 18 


proteins each (Table 5) 


Table 5 Gene set Enrichment analysis of modulated inflammatory pathways 


disease 


ID Path puaeealte 
athways ene code 
iis Gene Set 
: CASP7, NDUFA4, CASP3, PRKAB1, CXCL8, 
Non-alcoholic fatty liver 
hsa04932 15 RXRA, CASP8, GSK3B, IL6, AKT1, TNF, PPARA, 


NFKB1, ERN1, NR1H3 
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hsa05310 | Asthma 1 TNF 
hea0A064 NF -kappa B signaling 9 PTGS2, CXCL8, NFKBIA, BIRC2, TNF, MALT1, 
sa 
pathway NFKB1, LCK, CSNK2A1 
hsa05321 | Inflammatory bowel disease | 6 NFKB1, TNF, RORC, RORA, IL2, IL6 
hsa05323 | Rheumatoid arthritis 7 CTSK, CTSL, CXCL8, IL6, TNF, VEGFA, MMP1 
Arrhythmogenic right 
hsa05412 : . 4 ITGA2B, ITGB3, CACNAIC, ATP2A1 
ventricular cardiomyopathy 
Chemokine signaling CXCL8, NFKBIA, GSK3B, AKT1, CXCR1, NFKB1, 
hsa04062 7 
pathway PIK3CG 
hsa04145 | Phagosome 3 MPO, ITGB3, CTSS, CTSL 
hsa04514 | Cell adhesion molecules 1 PTPRC 
. ia IL6, NFKBIA, RORC, RXRA, RORA, AHR, IL2, 
hsa04659 | Th17 cell differentiation 10 
NFKB1, HSP9OAA1, LCK 
: . CASP7, CASP3, PTGS2, NFKBIA, BIRC2, CASP8, 
hsa04668 | TNF signaling pathway 11 
TNF, MMP9, NFKB1, IL6, AKT1 
Leukocyte transendothelial 
hsa04670 ; ; 3 PRKCA, MMP2, MMP9 
migration 
RPS6KA3, PRKAB1, NFKBIA, NOS3, GSK3B, IL6, 
hsa04931 | Insulin resistance 13 AKT1, TNF, PPARA, NFKB1, PTPN1, NR1H2, 
NR1H3 
Systemic lupus 
hsa05322 4 GRIN2B, ELANE, TNF, C5 
erythematosus 
hsa05416 | Viral myocarditis 5 ABL1, ABL2, CASP3, CASP8, CASP9 
hsa04614 | Renin-angiotensin system 3 CMA, REN, ACE 
AGE-RAGE signaling 
. NOS3, CASP3, PIM1, PRKCA, CXCL8, IL6, AKT1, 
hsa04933 | pathway in diabetic 12 
ane TNF, VEGFA, NOX4, NFKB1, MMP2 
complications 
PGF, IGF1R, RPS6KA3, CASP3, AKT1, PRKCA, 
hsa04010 | MAPK signaling pathway 18 CACNAIC, EGFR, NTRK2, FLT3, MET, TNF, 
VEGFA, NFKB1, KIT, PTPN7, FGF2, FGF1 
healS418 Fluid shear stress and is ITGB3, NOS3, CTSL, ITGA2B, NFKB1, AKT1, 
sal 
atherosclerosis GSTP1, TNF, VEGFA, MMP9, HSP90DAA1, MMP2 
. . IGF1R, NOS3, PRKCA, EGFR, NOS2, IL6, AKT1, 
hsa04066 | HIF-1 signaling pathway 9 
VEGFA, NFKB1 
Complement and 
hsa04610 : 4 F2, F10, C5, F7 
coagulation cascades 
: ; PTGS2 , CASP3, NFKBIA, CASP8, TNF, MMP9, 
hsa04657 | IL-17 signaling pathway 12 
NFKB1, CXCL8, GSK3B, MMP1, IL6, HSP90AA1 
Th1 and Th2 cell NFKB1, NFKBIA, IL2, LCK 
hsa04658 . pac 4 
differentiation 
Fc epsilon RI signaling 
hsa04664 4 PRKCA, AKT1, ALOX5, TNF 
pathway 
CASP7, NDUFA4, ATP2A1, BACE1, PSENI1, 
: . CACNAIC, CASP8, CASP9, GRIN2B, TNF, 
hsa05010 | Alzheimer disease 17 
CAPN2, CAPN1, CDK5, GRIN1, GSK3B, ERN1, 
CASP3 
Staphylococcus aureus 
hsa05150 | . : 1 C5 
infection 
hsa05163 | Human cytomegalovirus 18 PTGS2, ITGB3, CASP3, PTGER4, CXCL8, PRKCA, 
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infection NFKBIA, EGFR, CASP8, PTGER3, IL6, AKT1, 
TNF, VEGFA, GSK3B, CDK6, NFKB1, CASP9 
Human T-cell leukemia NFKBIA, IL6, AKT1, CDK2, TNF, IL2, NFKB1, 
hsa05166 ; . ; 9 
virus 1 infection TERT, LCK 
: . PTGS2, CASP3, CXCL8, NFKBIA, CASP8, CASP9, 
Kaposi sarcoma-associated 
hsa05167 : . . 14 IL6, AKT1, VEGFA, GSK3B, CDK6, NFKB1, 
herpesvirus infection 
PIK3CG, FGF2 
HDAC1, HDAC3, HDAC2, HDAC5, CASP3, 
hsa05203 | Viral carcinogenesis 16 HDAC6, HDACI, HDAC8, NFKBIA, CASP8, 


CDK1, CDK2, CDK6, NFKB1, HDAC10, HDAC4 
HDAC1, PTGS2, ITGB3, CYP24A1, CYP1B1, 
hsa05206 | MicroRNAs in cancer 17 CASP3, PIM1, PRKCA, EGFR, MET, VEGFA, 
CDkK6, MMP9, ABL1, NFKB1, ABCC1, HDAC4 
IGF1R, AKT1, PRKCA, EGFR, GSK3B, MET, 
GSTP1, CDK6, TERT, TXNRD1 


hsa05225 | Hepatocellular carcinoma 10 


Choline metabolism in 
hsa05231 3 PRKCA, AKT1, EGFR 
cancer 


Protein-protein network 

The protein-protein interaction (PPI) network is shown in Figure 3. RAC-alpha serine/threonine protein kinase (AKT1) and heat 
shock protein HSP 90-alpha (HSP90AA1) had the highest interactions, interacting with 31 proteins each at a 70% confidence level. 
The levels of AKT1, HSP90AA1 and TNF-a were e express. 


Network pharmacology analysis 

Seven hundred and ninety-seven (797) edges were identified in the component-target-pathway network in which five hundred and 
fourteen (514) were target-components interactions and two hundred and eighty-three (283) were target-pathway interactions 
(Figure 4). The constructed network included one hundred and ninety-six (197) nodes representing sixty-seven (67) 
phytoconstituents, ninety-eight (97) targets and thirty-three (33) pathways. Prothrobin scored the highest edge count; interaction 
was found with thirty-one (31) protein molecules. Similarly, the MAPK signaling pathway and human cytomegalovirus infection 
modulated the highest number of proteins (eighteen proteins apiece). Rutin, beta-sitosterone, astragalin, beta-amyrin, quercetin and 
4-hydroxymellein are the most important nodes of the entire network. 


Molecular docking analysis 

The phytoconstituents of N. sativa and M. oleifera interacted with the selected inflammatory targets at different levels of binding 
affinity. Table 6 shows the binding affinities of different compounds selected based on the results from the network analysis and the 
different proteins selected. The binding energies of most of the compounds were lower than -5Kcal/mol. To further validate the 
results of network pharmacology, most of the phytoconstituents of N. sativa and M. oleifera were docked with the predicted targets 
in addition to interferon-gamma. The results showed that some of the compounds had lower binding affinities than the ones 
predicted by network pharmacology. However, most of the compounds predicted by network pharmacology interacted with more 
amino acid residues and formed more hydrogen bonds with the amino acid residues. Table 7 shows the binding affinities of most of 
the compounds from N. sativa and M. oleifera. Figures 5-15 show the three-dimensional (3D) and two-dimensional (2D) 
representations of the molecular interactions of the amino acid residues of the targets with the compounds predicted to have the 
lowest binding energies, with Figures 5-14A representing those predicted by network pharmacology and figures 5-14B representing 
those predicted after docking with most of the compounds. 
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Figure 3 Protein-protein interaction of regulated proteins by the phytoconstituents from N. sativa and M. oleifera. 
Each node (circle) represents a protein and each line refers to an interaction. The thickness of the line reflects the strength of the protein-protein interaction. Circles without lines linking them to other 
circles indicate that the protein interacted with no other protein. 
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Figure 4 Components-targets—pathways network of phytoconstituents of N. sativa and M. olifera 


Green, pink and blue colours indicate the phytoconstituents, targets and pathways respectively. The lines (edges) represent the interactions. 


Table 6 Binding affinities of compounds of N. sativa and M. oleifera against the selected targets 


AG Energy (Kcal/mol) 
. PubChem 

Phytoconstituents HSP90 PTGS 

CID NR1H3 | F2 AKTI1 ree IL2 | NFKB1 5 ALOX5 | TNF IL6 | IFN-y 
24-methylene 

94204 -11.2 -8.4 -5.7 -8.1 -5.3 | -6.4 7A -8.1 -6.3 -6.4 | -8.9 
cycloartenol 
24-methylene- 

92113 -8.9 -7.5 -6.7 -8.6 -5.5 | -6.7 -7.6 -8.2 -6.4 -6.1 | -8.3 
cholesterol 
28-isoavenasterol 

91746804 -8.1 7.5 -5.6 -7.2 -5.3 | -5.9 -7.2 -8.1 -5.2 -6.3 | -8.4 
acetate 
4-hydroxymellein 169539 -7.6 -5.7 -4.8 -6.7 45 | -5.7 -6.8 -74 -4,9 -5.5 | -5.8 
5-dehydro- 

44263331 -8.9 -8.2 -6.6 -8.7 -5.8 | -7.1 -7.9 -9.0 -6.6 -5.9 | -8.8 
avenasterol 
Alpha-spinasterol 5281331 -8.8 -8.3 -6.1 -8.6 -5.8 | -6.8 -9.2 -9.6 -6.2 -7.0 | -8.9 
Astragalin 5282102 -9.1 -7.8 -5.4 -8.8 -5.4 | -7.0 -8.4 -7.9 -5.7 -6.4 | -6.8 
Benzyl-amine 7504. -5.8 -4.6 -3.7 -4.9 -3.1 | -43 -5.2 -5.2 -3.8 3.5 | -4.3 
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Beta-amyrin 73145 -9.7 -8.7 -6.5 -7.6 -6.5 | -7.8 -9.2 -10.0 -6.8 -6.5 | -10.7 
Beta-sitosterone 9801811 -8.1 -7.6 -6.7 -8.4 -5.8 | -6.1 -7.3 -8.6 -6.1 -6.2 | -8.0 
Brassicasterol 5281327 -12.0 7.9 -6.7 -7.7 -5.5 | -6.6 -9.0 -9.2 -6.7 -6.1 | -8.6 
Butulic acid 64971 9.2 7.4 -5.9 -6.6 -5.8 | -6.7 -8.1 -9.4 -6.6 -6.3 | -8.4 
Butyrospermol 12302182 -10.0 -8.6 -6.0 -8.8 -5.7 | -6.9 -7.5 -9.1 -5.5 -6.9 | -9.8 
Caffeic acid 689043 7.1 -6.1 4.2 ~5.9 -4.7 | -5.1 -6.9 -6.8 -4.7 -5.1 | -6.1 
Campestanol 119394 -11.2 -7.3 -6.4 -8.3 -5.3 | -6.5 -7.5 -8.4 -6.4 -6.1 | -8.3 
Campesterol 173183 -10.0 -8.0 -6.5 -8.5 -5.4 | -6.6 -8.8 -8.4 -6.5 -5.8 | -9.1 
Carcacol 10364 ~7.2 -5.5 4.3 -5.9 -44 | -4.6 -6.6 -7.0 -4.2 -44 | -5.7 
Carvone 7439 -6.9 -5.2 4.1 -5.8 4.2 | -4.5 -6.6 -6.5 -4.4 -4.3 | -5.4 
Cholesterol 5997 -11.6 -7.8 -6.3 -8.1 5.4 | -6.4 -8.6 -79 -6.0 5.8 | -8.2 
Citrostadienol 9548595 -11.6 -8.0 -5.8 -8.9 -5.8 | -6.3 -9.3 -8.5 7.1 -6.7 | -9.5 
Clerosterol 5283638 -11.5 7.4 -6.3 -8.3 -5.6 | -6.2 -7.2 -8.9 -6.0 -5.8 | -8.2 
Cycloeucalenol 101690 -8.2 -8.2 -5.6 -8.3 -5.9 | -6.8 -9.2 -8.6 -6.5 -6.7 | -8.6 
Damascenine 21368 -6.7 -5.2 3.9 -5.5 -3.8 | -4.6 -6.3 -6.3 -4.2 -4.3 | -5.0 
Dehydroascorbic acid | 440667 -5.5 -5.8 4.3 -6.1 -4.1 | -5.0 -6.5 -6.4 -49 49 | -4.6 
Delta-5-avenasterol 5281326 -9.6 7.5 -6.7 -8.1 5.4 | -6.4 -7.6 -8.2 -6.9 -6.0 | -9.6 
Delta-7-avenasterol 12795733 -11.0 -8.5 -6.1 -8.8 -5.5 | -7.0 -8.4 -8.9 -6.6 -6.6 | -8.9 
Delta-tocopherol 92094 8.9 -6.0 4.3 -7.0 4.5 | -4.9 -7.0 -7.8 -4.3 -5.3 | -7.7 
D-galactose 6036 -5.8 4.8 -4.4 -5.6 -3.8 | -4.7 -6.2 -6.1 -44 -4.7 | -5.1 
Dithymoquinone 398941 7.9 7.3 -5.3 -6.7 -5.1 | -6.3 -7.5 -8.4 -5.7 -6.0 | -7.4 
D-xylose 135191 -5.3 5.9 4.3 -5.0 3.5 | -4.5 -5.7 -5.3 -4.1 -3.9 | -4.6 
Gamma-tocopherol 92729 7.6 -6.3 --43 -6.8 -4.6 | -5.2 -7.6 -6.8 -4.8 -5.6 | -8.3 
Gramisterol 5283640 -8.0 7.4 -6.3 -8.5 5.5 | -6.3 -8.4 -8.6 -6.6 -6.6 | -8.7 
Hederagenin 73299 -8.8 -7.7 -6.1 -7.3 6.2 | -7.2 -8.4 -9.1 -6.2 -6.3 | -9.8 
Indole-3-acetic acid 802 7.4 -5.8 -4.5 -6.7 -4.5 | -4.9 -6.8 -7.0 -4.8 -5.0 | -5.7 
Indoleacetonitrile 351795 7.1 -5.7 4.3 -6.4 3.9 | -5.1 -6.4 -6.8 -4.6 -4.7 | -5.9 
Kaempferol 5280863 -8.9 7.4 -5.7 7A 5.5 | -6.4 -9.3 -8.0 -5.8 -6.0 | -6.7 
Lophenol 160482 -11.3 -7.7 -5.9 -8.2 -5.2 | -6.4 -7.0 -8.5 -6.1 -6.4 | -8.7 
L-rhamnose 25310 -5.4 4.9 4.3 -5.4 -3.7 | -4.6 5.9 -5.7 -4.0 -4.5 | -4.2 
Maltotriose 439586 7.2 -7.0 -4.8 -7.7 -5.1 | -6.3 -6.7 -7.6 -5.1 -5.4 | -5.6 
N-acetyl-L-cysteine 12035 -5.3 4.5 -3.4 -5.0 -3.4 | -4.0 4.9 -5.0 -3.9 -3.8 | -3.9 
Obtusifoliol 65252 -8.1 7.9 -5.9 -8.2 5.3 | -6.9 -8.1 -8.9 -6.3 -6.6 | -9.2 
Octadeca-9,12-dienoic | 3931 7.2 -5.1 4.4 -5.6 -3.3 | -4.7 -5.8 -5.8 -3.9 4.6 | -5.8 
acid 
Oleic acid 445639 -6.9 -5.3 -3.6 -5.2 -3.8 | -4.0 -5.6 -5.8 -3.9 -4.2 | -5.6 
Phytosterols 12303662 -8.3 -8.1 -6.4 -8.6 -5.3 | -6.1 ~7.2 -8.8 -6.2 -5.6 | -8.1 
Quercetin 5280343 -9.1 -8.1 -5.5 -7.5 -5.6 | -6.7 -9.4 -8.3 -5.7 -6.4 | -6.6 
Quercetin-3’- 25203368 -8.8 -8.5 -6.1 -8.5 -5.4 | -6.6 8.9 7.9 -6.1 -6.7 | -6.8 
glucoside 
Rutin 5280805 -9.3 -8.3 -6.5 -7.7 -6.2 | -7.5 -10.0 | -8.6 -6.8 -6.6 | -7.7 
Stigmast-7-en-3-beta- | 441837 -11.1 -8.2 -5.7 -7.7 -5.7 | -6.7 7.4 -8.4 -6.3 -6.7 | -8.6 
ol 
Stigmastanol 241572 -11.1 -7.6 -6.3 -8.4 -5.2 | -6.1 7.3 -8.3 -6.1 -6.1 | -8.0 
Taraxerol 92097 -10.0 -8.1 7.2 -7.3 -6.8 | -7.6 9.1 -9.9 -6.5 -6.7 | -10.1 
Thymohydroquinone | 95779 -6.7 -5.5 4.1 -6.2 -4.0 | -4.7 -6.3 -7.1 -4.3 -4.8 | -5.4 
Thymol 6989 -6.7 -5.4 4.1 -5.8 -4.0 | -4.8 -6.0 -7.1 -4.0 -4.7 | -5.7 
Thymoquinone 10281 -7.0 -5.5 -4.0 -6.4 -3.9 | -4.6 -6.7 -7.3 -4.3 -4.7 | -5.3 
Tirucallol 101257 -9.2 -7.2 7.1 -7.8 -5.4 | -6.7 -8.6 -9.5 -6.0 -7.0 | -9.5 
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Vanillin 1183 -6.4 -5.4 3.9 -5.6 -3.9 | -4.6 -6.1 -6.1 -4.2 4.1 | -4.5 
*Diclofenac 3033 7.9 -6.7 -5.2 -79 -4.7 | -5.7 7.7 -7.1 -5.0 -5.1 | -6.1 

Table 7 Binding affinities of the compounds against the selected targets based on network pharmacology 

PubChem | AG Energy (Kcal/mol) 
Compounds 
CID NR1H3 F2 AKT1 | HSP90AA1 | IL2 | NFKB1 | PTGS2 | ALOX5 | TNF | IL6 

24-methylene cycloartanol 94204 -113 Ss 5 - - - = 7 - : 
24-methylene-cholesterol 92113 -6.8 8.5 | - - - - 5 - - 
28-isoavenasterol acetate 91746804 | -6.9 = - - - - - - - = 
4-hydroxymellein 169539 - - : - - = * = -5.1 
5-dehydro-avenasterol 44263331 | -7.0 -8.0 | - - - - - - - 
Alpha-spinasterol 5281331 -7.0 -7.8 | - - - - . ; - 
Astragalin 5282102 - -8.2 | - = 5.5 | - - -8.0 -5.7 - 
Benzyl-amine 7504 5 45 | - ; - - : ss = ; 
Beta-amyrin 73145 -7.7 -8.7 | - 7 = -7.5 - 7 7 . 
Beta-sitosterone 9801811 -6.7 : - -7.1 - - = - 7 
Brassicasterol 5281327 -7.2 7.7 | - - - - - - - - 
Butyrospermol 12302182 | -10.6 7.3 | - = - -7.0 - - - - 
Campestanol 119394 “11.2 “ - - - - : - 
Campesterol 173183 -12.3 8.2 | - - - - - - - - 
Cholesterol 5997 -11.7 -7.9 | - = - - - - - - 
Citrostadienol 9548595 -7.6 -8.1 | - = - -6.6 - - “ ' 
Clerosterol 5283638 -6.8 7A | - - - - - - - - 
Cycloeucalenol 101690 -10.8 - - - - - - - - - 
Damascenine 21368 - - - - - - - -5.4 - - 
Dehydroascorbic acid 440667 - = - - - - “0.9 - 
Delta-5-avenasterol 5281326 -13.0 7A - 
Delta-7-avenasterol 12795733 | -11.3 8.5 | - - - - = : : 7 
Delta-tocopherol 92094 = = 4.3 = - - 7.7 ~7.2 - = 
D-galactose 6036 - 4.9 | - - - - - - - - 
D-xylose 135191 = -6.0 | - = - - - - - - 
Gamma-tocopherol 92729 -8.1 “ -4.5 - - “70 -6.3 - - 
Gramisterol 5283640 -11.5 7.9 | - - - -6.3 - “ - - 
Hederagenin 73299 “77 - - . - “7.1 -8.7 - - - 
Indole-3-acetic acid 802 - - - - - - -6.7 : - - 
Kaempferol 5280863 - 7.2 | - - - - 7.7 -8.2 - = 
Lophenol 160482 -11.7 7.6 | - - - -6.4 - - - = 
L-rhamnose 25310 = 49 | - - - - . . = 
Maltotriose 439586 - 7.0 | - - - - - - - - 
Obtusifoliol 65252 -10.7 7.8 | - - - - - - - - 
Phytosterols 12303662 | -6.7 -8.0 | - - - - - - - - 
Quercetin 5280343 - 7.2 | - = - - -9.5 -8.3 - = 
Quercetin-3’-glucoside 25203368 - -7.6 | - - 5.4 | - - -7.6 -5.8 - 
Rutin 5280805 = -8.7 | - - “6.2 | - - -8.6 -6.9 - 
Stigmast-7-en-3-beta-ol 441837 -7.2 -8.2 | - - - - - . “ - 
Stigmastanol 241572 -6.9 - - - - - - “ - 
Taraxerol 92097 - -8.1 | - - - 7.3 - - - 
Thymol 6989 - - - - - - - 7.1 - - 
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Thymoquinone 10281 - - - - - - -6.1 - - 
Tirucallol 101257 -75 -7.8 - - - - - - - 
*Diclofenac 3033 - - - - - -7.7 - - - 
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Figure 5 Molecular interactions of amino-acid residues of NR1H3 with A. Delta-5-avenasterol and B. Brassicasterol showing the 2D 


(right) and 3D (left) views 


A. Delta-5-avenasterol interacted with different amino-acid residues in the binding site of NR1H3 forming pi-sigma, alkyl, pi-alkyl bonds as well as van der Waals forces. B. Brassicasterol reacts with 


different amino acid residues at the binding site of NR1H3 forming alkyl and pi-alkyl bonds and van der Waals forces with the protein. 
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Figure 6 Molecular interactions of amino-acid residues of F2 with A. Rutin and B. Beta-amyrin showing the 2D (right) and 3D (left) 


174 


B van der Waals O ay 
B Pi-Sigma CL] Pi-Alkyl 


views 
A. Rutin interacted with different amino acid residues present in the binding site of F2 forming different kinds of bonds. Rutin formed hydrogen bonds with ARG35, GLU39, LEU40, TRP60D, ARG73 
and ASN143. B. Beta-amyrin reacts with different amino acid residues at the binding site of F2 forming alkyl, pi-alkyl, pi-sigma bonds and van der Waals forces with the protein. 
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Figure 7 Molecular interactions of amino-acid residues of AKT1 with A. Delta-tocopherol and B. Taraxerol showing the 2D (right) 


and 3D (left) views 
A. Delta-tocopherol interacted with VAL7, ASN31, ASP32, THR105, VAL106, GLY109, LEU110, LYS112 and GLN113 present in the binding site of AKT1 forming different kinds of bonds with the 


protein. B. Taraxerol reacts with different amino acid residues at the binding site of F2 forming alkyl bonds and van der Waals forces with the protein. 
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Figure 8 Molecular interactions of amino-acid residues of HSP90AA1 with A. Beta-sitosterone and B. Citrostadienol showing the 2D 
(right) and 3D (left) views 


A. Beta-sitosterone interacted with different amino acid residues present in the binding site of HSP90AA1 forming alkyl bonds and van der Waals forces with the protein. B. Citrostadienol reacts 
with different amino acid residues at the binding site to HSP90AA1 forming alkyl and pi-alkyl bonds and van der Waals forces with the protein. alkyl bonds and van der Waals forces with the 


protein. 
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Figure 9 Molecular interactions of amino-acid residues of IL2 with A. Astragalin and B. Taraxerol showing the 2D (right) and 3D 
(left) views 


der Waals forces (GLU67, ASN71, GLN74, ASN77 and HIS79). 


A. Astragalin formed two hydrogen bonds at the target site of IL2, one of which is a conventional hydrogen bond with ASN71. Astragalin also formed a pi-alkyl bond with ARG81 and van der 
Waals forces with different amino acid residues. B. Taraxerol binds at the target site of IL-2 with one conventional hydrogen bond (PHE78), alkyl bonds with ARG81, PRO82 and ARG83 and five van 
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Figure 10 Molecular interactions of amino-acid residues of NFKB1 with beta-amyrin showing the 2D (right) and 3D (left) views 
Beta-amyrin reacts with different amino acid residues at the binding site of NFKB1 forming alkyl bonds and van der Waals forces 
with the protein. 
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Figure 11 Molecular interactions of amino-acid residues of PTGS2 with A. Quercetin, B. Rutin and C. Diclofenac (standard drug) 


showing the 2D (right) and 3D (left) views 

A. Quercetin binds firmly at the target site of PTGS2 with four conventional hydrogen bonds (CYS36, GLY45, TYR130 and GLU465), one carbon-hydrogen bond (VAL46), one pi-sigma bond 
(LEU152), two pi-alkyl bonds (CYS47 and PRO153) and twelve van der Waals forces. B. Rutin binds firmly at the target site of PTGS2 with eight conventional hydrogen bonds (CYS47, ASP58, 
THR60, GLN61, ASP125, TYR130, LYS137), three carbon-hydrogen bonds, one pi-sigma bond, three pi-alkyl bonds and ten van der Waals forces C. Diclofenac bind firmly at the target site of PTGS2 
with two conventional hydrogen bonds (ARG44 and GLY45), six alkyl and pi-alkyl bonds and seven van der Waals forces. 
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Figure 12 Molecular interactions of amino-acid residues of ALOX5 with A. Rutin and B. Beta-amyrin showing the 2D (right) and 3D 


(left) views 

A. Rutin binds firmly at the target site of ALOX5 with four conventional hydrogen bonds (ARG101, GLN129, LYS133 and GLU134), one carbon-hydrogen bond (HIS130), one pi-pi T-shaped alkyl 
interaction (HIS130), pi-alkyl bonds with VAL110 and LYS133 and five van der Waals forces. B. Beta-amyrin binds at the target site of ALOX5 with three alkyl bonds (ARG143, TRP144 and ILE330), 
one pi-sigma bond (TYR515) and eleven van der Waals forces. 
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Figure 13 Molecular interactions of amino-acid residues of TNF-a with A. Rutin and Citrostadienol showing the 2D (up) and 3D 


(down) views 
A. Ruin binds firmly at the target site of TNF-a with three conventional hydrogen bonds (PRO100, GLN102 and LYS112), pi-pi T-shaped bonds with TRP114, one pi-alkyl bond (CYS69) and nine van 
der Waals forces. Citrostadienol binds at the target site of TNF-a forming alkyl and alkyl interactions with five amino acids (PRO70, PRO106, ALA111, TYR141 and PHE144) and van der Waals 


forces with nine amino acids. 
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Figure 14 Molecular interactions of amino-acid residues of IL-6 with 4-hydroxymellein showing the 2D (right) and 3D (left) views. 


A. 4-hydroxymellein binds firmly at the target site of IL- with two conventional hydrogen bonds (GLU95 and LYS120), one pi-anion 
interaction with GLU95, one alkyl bond and six van der Waals forces. B. Tirucallol binds at the target site of IL-6 forming a 
hydrogen bond with GLU99, alkyl and alkyl interactions with five amino acids and van der Waals forces with six amino acids. 
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Figure 15 Molecular interactions of amino-acid residues of IFN-y with beta-amyrin showing the 2D (right) and 3D (left) views 
Beta-amyrin binds at the target site of IFN-y with ten alkyl and pi-alkyl bonds and two van der Waals forces (PHE57 and LYS80). 


| van der Waals | Pi-Akyl 


4. DISCUSSION 


Uncontrolled inflammation could eventually result in disorders such as rheumatoid arthritis, psoriasis, periodontitis, inflammatory 
bowel disease and even cancers (Diakos et al., 2014). Non-steroidal anti-inflammatory drugs (NSAIDs) and corticosteroids are 
widely applied therapeutic agents which have achieved good therapeutic effects. Some side effects have however resulted the to 
their extensive use. Meanwhile, based on continuous breakthroughs in studies for compatibility of traditional medicine, increasing 
attention has been paid because of its significant therapeutic effects and minimal side effects (Li et al., 2019). 

Several in vivo and in vitro studies on Nigella sativa and Moringa oleifera reveal that these plants possess anti-inflammatory 
activity (Chehl et al., 2009; Tshabalala et al., 2019). Although research on the anti-inflammatory potential of these plants continues, 
their molecular mechanism of action remains unclear. Fortunately, network pharmacology, which is especially suitable for multi- 
compound and multi-target analysis, provides a prospective method to solve this problem with bioinformatics, molecular biology 
and databases (Jiang et al., 2021). 

In this study, network pharmacology and molecular docking were used to analyze the anti-inflammatory mechanism of Nigella 
sativa and Moringa oleifera. Network analysis predicted the modulation of ninety-eight inflammatory proteins by sixty-four 
phytoconstituents which suggests multiple compounds and protein interaction. Based on network pharmacology data, 
phytoconstituents, including rutin, beta-sitosterone, astragalin, beta-amyrin, quercetin and 4-hydroxymellein and other compounds 
are the most important nodes of the entire network. Previous studies have shown that rutin, beta-sitosterone, astragalin, beta- 
amyrin, quercetin and others, have an obvious anti-inflammatory effect (Yoo et al., 2014; Paniagua-Perez et al., 2017; Hu et al., 2020; 
Melo et al., 2011; Saeedi-Boroujeni and Mahmoudian-Sani, 2021). PTGS2, TNF-a, IL6 and NFKBlare the targets associated 
connected with inflammation-related pathways such as TNF and IL-17 signaling pathways. These results suggest that N. sativa and 
M. oleifera could probably regulate inflammation as well as the immune system. PTGS2, also known as COX-2 is an enzyme that 
catalyzes arachidonic acid to produce prostaglandins, thereby playing a pro-inflammatory role and unsettling the balance of the 
internal environment (Jain et al., 2015). 

PTGS2 can be inhibited by inhibitors, thus blocking the conversion of arachidonic acid into prostaglandins, thus reducing the 
level of prostaglandins and alleviating inflammation and pain symptoms of tissues involved in prostaglandins (Zhang et al., 2011). 
TNF-a and IL-6 are pro-inflammatory cytokines that are involved in the upregulation of inflammatory reactions (Zhang and An, 
2007). TNF-a can mediate inflammatory reactions, induce the production of IL-18, IL-6 and other inflammatory factors and 
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aggravate the local inflammatory reaction (Turner et al., 2007). IL-6 is produced locally in the inflamed tissue following cellular 
activation by bacterial lipopolysaccharide or other cytokines such as IL-16 or TNF-a (Roshene and Ramesh, 2017). TNF-a is mainly 
secreted by activated T lymphocytes and mononuclear macrophages. On the other hand, IL6 is a signaling molecule that can be 
secreted by various cells of the body including, lymphocytes, macrophages/monocytes and epithelial cells (Zhu et al., 2021). NFKB1 
(nuclear factor — kappa B1) also known as p50 is an inducible transcription factor that regulates a large number of genes involved in 
diverse processes of the immune and inflammatory responses (Liu et al., 2017). 

The target levels of AKT1, HSP90AA1 and TNF-a were excessively expressed in the protein-protein interaction (PPI) network 
suggesting that they may play a pivotal role in the anti-inflammatory effects of N. sativa and M. oleifera. AKT1 (RAC-alpha 
serine/threonine-protein kinase) is a member of the AKT kinase family that regulates metabolism, proliferation, cell survival, gene 
expression, growth and angiogenesis through a series of downstream substrates (Tianyu and Liying, 2021). Studies reveal that 
AKT1 is essential for acute inflammation and exerts its effects primarily through the regulation of vascular permeability, leading to 
edema and leukocyte extravasation (Di Lorenzo et al., 2009). The heat shock protein 90a (HSP90a), encoded by the HSP90AA1 gene 
is a 90kDa ATP-dependent molecular chaperone that guides the late-stage tertiary folding and maintains the conformational 
integrity of multiple proteins such as oncogenic proteins, steroid receptors and transcription factors (Taipale et al., 2010). 

HSP90 is widely expressed in eukaryotic cells but usually in a latent, uncomplexed form whereas tumors express high levels of 
catalytically active HSP90 found in complex with oncogenic client proteins. HSP90 is also highly expressed in inflammatory cells 
(Lilja et al., 2021). HSP90 inhibitors have been shown to ameliorate inflammatory and auto-immune diseases such as 
encephalomyelitis, rheumatoid arthritis and systemic lupus erythematosus (Tukaj and Wegrzyn, 2016). Similarly, thirty-three 
pathways were identified to be associated with proteins related to inflammation with reference to KEGG pathway analysis. The key 
targets were mainly concentrated in pathways such as the TNF signaling pathway, Human cytomegalovirus infection, and MAPK 
signaling pathway. The main proteins involved in these pathways include TNF-a, PTGS2, IL2, IL6é, NFKB1, CASP3, and NFKBIA 
among others. 

Furthermore, molecular docking technology was used to verify the core active components and target proteins on the prediction 
results of network pharmacology. The proteins F2, NR1H3, AKT1, HSP90AA1, IL2, NFKB1, PTGS2, ALOX5, TNF-a and IL6 were 
docked with the different compounds predicted by the network pharmacology. The verification results showed that the binding 
free energies of the compounds with the targets were much less than -5 Kcal/mol, which indicated that the ligand molecules could 
bind spontaneously (Abdullahi and Adeniji, 2020). To further validate the docking results, molecular docking was carried out 
between most of the compounds of N. sativa and M. oleifera and the aforementioned proteins in addition to IFN-y (Interferon- 
gamma). The results showed that some of the compounds had lower binding affinities than the ones predicted by network 
pharmacology. 

However, the compounds predicted by network pharmacology (apart from the interaction between rutin and PTGS2) interacted 
with more amino acid residues and formed more hydrogen bonds with the amino acid residues. Hydrogen bonds are important in 
protein-ligand interaction because they stabilize the ligand in the binding pocket (Kostal, 2016). Hydrogen bonds occur between a 
hydrogen atom’s negatively charged groups (e.g., N, O, F). When several hydrogen bonds occur simultaneously at the binding site, 
they increase the strength and stability of drug-receptor interaction substantially. However, too many hydrogen bond 
donors/acceptors (like in the case of rutin and PTGS2 interaction) can have a deleterious effect on the drug’s membrane partition 
and permeability. Too many hydrogen bonds can reduce the affinity toward the hydrophobic membrane region and increase the 
desolvation penalty upon penetration of the drug (Coimbra et al., 2021). All of this gives credence to the fact that the network 
pharmacology tool is effective in drug development. 

IFN-y also had a strong binding affinity with some of the compounds, including beta-amyrin (-10.7Kcal/mol), taraxerol (- 
10.1Kcal/mol) and hederagenin (-9.8Kcal/mol) among others. Drug-likeness character of each phytoconstituents from N. saiva and 
M. oleifera was predicted based on “Lipinski’s Rule of Five” which is a qualitative assessment to evaluate the drug-like property 
concerning factors like oral bioavailability (Lipinski, 2004). In the present study, most of the compounds possess positive drug 
likeness scores reflecting their oral bioavailability. Similarly, the ADMET of each phytoconstituents was predicted to assess the 
important pharmacokinetic parameters in the human body including the probable toxicity of each molecule and was compared 
with diclofenac. 


5. CONCLUSION 
In conclusion, a network-based pharmacological analysis was utilized to predict the molecular targets for the anti-inflammatory 
activity of Nigella sativa and Moringa oleifera. Target proteins, including NR1H3, F2, AKT1, HSP90AA1, IL2, NFKB1, PTGS2, ALOX5, 
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TNF-a, IL6 and IFN-y, as well as signaling pathways such as TNF, MAPK, IL-17 and HIF-1 were linked to the anti-inflammatory 
activity of N. sativa and Moringa oleifera. Based on network pharmacology data and molecular docking, phytoconstituents which 


include rutin, beta-sitosterone, astragalin, beta-amyrin and quercetin were predicted to inhibit key inflammatory proteins. These 


molecular targets may contribute to the overall anti-inflammatory activity of N. sativa and M. oleifera and its potential therapeutic 


applications for several inflammatory-mediated diseases. 
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